Thermal effects on quark-gluon mixed condensate g(qa fJjU G^ u q) from lattice QCD 
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We present the first study of the thermal effects on the quark-gluon mixed condensate 
giqo-^vG^vq) , which is another chiral order parameter, in SU(3) C lattice QCD with the Kogut- 
Susskind fermion at the quenched level. Using the lattices at (3 = 6.0, 6.1 and 6.2 in high statistics, 
we calculate g{qc ' ^q) as well as (qq) in the chiral limit for & T 500MeV. Except for the 
sharp decrease of both the condensates around T c ~ 280MeV, the thermal effects are found to be 
weak below T c . We also find that the ratio m\ = g(qa ^„G M „q) / '{qq) is almost independent of the 
temperature even in the very vicinity of T c . This result indicates nontrivial similarity in the chiral 
behaviors of the two different condensates. 

PACS numbers: 12.38.Gc, 12.38.-t, 11.15.Ha 



I. INTRODUCTION 

Quantum chromodynamics (QCD) exhibits interesting 
nonperturbative phenomena such as spontaneous chiral- 
symmctry breaking and color confinement, which are 
considered to be related to the nontrivial structure of 
the QCD vacuum. In order to clarify the mechanism of 
these phenomena and their relation to the QCD vacuum, 
extensive studies have been made. In fact, at high tem- 
perature, QCD is believed to exhibit phase transition into 
the quark-gluon plasma (QGP), where chiral symmetry 
is restored and the color is deconfined. The QGP phase 
is considered to have been actually realized in the early 
universe, and the on-going RHIC experiments attempt 
to produce QGP in the laboratory through relativistic 
heavy- ion collisions, which motivates further studies of 
finite temperature QCD. 

For the theoretical study of the QCD vacuum struc- 
ture at finite temperature, the condensates such as (qq), 
a s (G M „G pi,) and g(qa ^ V G ^q) play the relevant role, be- 
cause they characterize the nontrivial QCD vacuum di- 
rectly. In fact, at finite temperature, changes of the struc- 
ture of the QCD vacuum will be represented by the ther- 
mal effects on the condensates. 

In this paper, we study the thermal effects on 
the quark-gluon mixed condensate g(qa M „G ^q) = 
g(qv \ivG^ v \\ A q) for the following reasons. First of all, 
we note that the chirality of the quark in the operator 
g^a^G^q) flips as 



g(qo ^vG ^ v q) 
= g(qR(o-^G^)q L ) + g^L^^G^qa). 



(1) 



Therefore, g(qa ^vG ^q) plays the role of a chiral order 
parameter, which can indicate chiral restoration at finite 
temperature, as well as the usual quark condensate (qq). 
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Secondly, we emphasize that g(qa p, v G ^q) character- 
izes different aspect of the QCD vacuum from (qq). In 
particular, g (qa M „G ^ v q) reflects the color-octet compo- 
nents of quark-antiquark pairs in the QCD vacuum, while 
(qq) reflects only the color-singlet q-q components. The 
mixed condensate g(qa M „G ^q) thus represents the direct 
correlation between color-octet q-q pairs and the gluon 
field strength G A V , i.e., the color-electromagnetic field 
spontaneously generated in the QCD vacuum [l| . There- 
fore, the thermal effects on g(qa^G^ u q) in comparison 
with that of (qq) will give new and important informa- 
tion on chiral restoration of the QCD vacuum at finite 
temperature. 

As the third point, we show how the mixed condensate 
affects the properties of hadrons. To clarify the physical 
meaning of the condensates, the QCD sum rule is a use- 
ful framework, in which the hadronic properties can be 
directly connected to the various condensates with the 
help of the dispersion relation Q. One of the exam- 
ples where the mixed condensate plays an important role 
is the QCD sum rule for baryons 0, 0- I n the actual 
calculation of such sum rules, especially in the decuplet 
baryons, the contribution from g(qa ^ V G ^ v q) amounts to 
the same magnitude as the leading contribution from 
(qq), and the mixed condensate has large effects on the 
N-A splitting |^. In terms of the chiral properties of 
baryons, the parity splittings of baryons stem from the 
chiral-odd condensates such as (qq) and g(qo~ ^ v G^q) 0. 
Furthermore, recent study [|| shows that the parameter 
g(so ilvG ^ u s) j (ss) is a key quantity for the prediction on 
the parity of the recently discovered penta-quark baryon, 
6+(1540) 0. The mixed condensate is also important 
in other sum rules such as fight-heavy meson systems Q , 
and exotic meson systems .9]. In these respects, the eval- 
uation of the thermal effects on g(qo-^ v G^ v q) is expected 
to give a useful input in QCD sum rules to investigate 
hadron phenomenology at finite temperature. 

For the analysis of the thermal effects on the mixed 
condensate g(qa ^ V G ^q) , we use lattice QCD Monte 
Carlo simulation, which is the direct and nonperturba- 
tive calculation from QCD. So far, the mixed condensate 
at zero temperature has been analyzed phenomenologi- 
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cally in the QCD sum rules [ljj. ^ n the lattice QCD, 
a pioneering work ^l| was done long time ago, but 
the result was rather preliminary because the simula- 
tion was done with insufficient statistics using a small 
and coarse lattice. Recently, new lattice calculations 
have been developed by our group 0] using the Kogut- 
Susskind (KS) fermion, and by another group using 
the Domain- Wall fermion. At finite temperature, how- 
ever, there has been no result on g(qo^ v G^ v q) except for 
our early reports |l4j| . Therefore, we present in this pa- 
per the first and intensive results of the thermal effects 
on giqa^G^q) as well as (qq), including the analysis 
near the critical temperature. 

This paper is organized as follows. In Sec. |nj we 
explain the formalism to evaluate the condensates. In 
Sec. IIIII we present the lattice QCD data, and discuss 
the physical implication of the results. Sec. II VI is devoted 
to the summary of the paper. 



II. FORMALISM 

In the calculation of g (qo n V G ^q) and {qq), we use the 
SU(3) C lattice QCD with the KS-fermion at the quenched 
level. We note that the KS-fermion preserves the explicit 
chiral symmetry for the quark mass m = 0, which is a 
desirable feature to study both the condensates of chiral 
order parameters. In order to evaluate the condensates, 
we calculate the SU(4) t flavor-averaged condensates on 
the lattice as 



and 



a 3 (qq) 



a 5 g(qa f _ ll/ G f _ l „q) 



(2) 



£ Tr[{qf(x)qf(x))a^G^ 



(3) 



f, 



Here, "Tr" refers to the trace over the spinor and the 
color indices, and {q* \y)q j \x)) denotes the Euclidean 
quark propagator of the /-th flavor. For the gluon field 
strength Gjf*, we adopt the clover- type definition on the 
lattice to eliminate 0(a) discretization error, 



G 



A lat 



-Tr 



h.c. 



where denotes the plaquette operator, and A (A = 
1,2,- •• ,8) is the color SU(3) Gell-Mann matrix. In 
the KS-fermion formalism, SU(4) / quark-spinor fields, q 
and q, are converted into single-component Grassmann 
KS-fields, x an d X, respectively, with the proper inser- 
tion of gauge- link variables which ensure the gauge co- 
variance [IJ. Note that the contraction of flavor and 
spinor indices corresponds to summation over the KS- 
hypercube |12| and suppresses 0(a) discretization error. 



The detailed formulations and the diagrammatic repre- 
sentations for the calculation of the condensates are given 
in Ref. E^. 

We perform Monte Carlo simulations with the stan- 
dard Wilson gauge action for fi = 2N c /g 2 — 6.0,6.1 
and 6.2. The lattice units are obtained as a -1 = 1.9, 
2.3 and 2.7GeV for (3 = 6.0, 6.1 and 6.2, respectively, 
where we use the lattice data [l5j with the string tension 
= 427MeV. 

In order to perform the calculation at various temper- 
atures, we use the following lattices, 

i) /3 = 6.0, V = 16 3 x N t (N t = 16, 12,10,8,6,4), 

ii) = 6.1, V = 20 3 x N t (N t = 20, 12, 10, 8, 6), 

iii) (3 = 6.2, V = 24 3 x N t (N t = 24, 16, 12, 10, 8), 
which correspond to <, T < 500MeV with the spatial 
volume of L 3 ~ (1.6 — 1.8fm) 3 . The critical temperature 
T c is obtained from the Polyakov-loop susceptibility [l^ 
in terms of the confinement/deconfinement phase transi- 
tion. From the Polyakov-loop susceptibility on the above 
three lattices, we obtain T c j^d ~ 0.64, or T c ~ 280MeV, 
which is consistent with Ref. 

We generate 100 gauge configurations for each lattice, 
where we pick up each configuration for every 500 sweeps 
after 1000 sweeps for the thermalization. In the vicinity 
of the phase transition point, namely, 20 3 x 8 at (3 = 6.1 
and 24 3 X 10 at j3 = 6.2, the fluctuations of the conden- 
sates get larger. We hence generate 1000 gauge config- 
urations on these two lattices for the accurate estimate. 
For the lattices at T > T c , we only use the gauge config- 
urations which are continuously connected to the trivial 
vacuum 17^ = 1. 

To obtain the condensates, we calculate the Euclidean 
propagator by solving the matrix inverse equations. 
Here, we use the current-quark mass of ma = 0.0105, 
0.0184 and 0.0263 {fi = 6.0), ma = 0.00945, 0.0162 
and 0.0234 {/3 = 6.1) and ma = 0.00770, 0.0134 and 
0.0192 {(3 = 6.2), which correspond to the physical mass 
of m ~ 20, 35 and 50MeV, respectively. For the KS- 
fields, x an d Xj the anti-periodic condition is imposed 
at the boundary in both the temporal and spatial di- 
rections. In each configuration, we measure the con- 
densates on 16 different physical space-time points at 
(3 = 6.0 and 2 points at (3 = 6.1 and 6.2. These points 
are taken so as to be equally spaced in the lattice 4- 
dimensional volume, i.e., on the lattices with the volume 
V = L 3 x N t = (2l) 3 x (2n t ) in the lattice unit, we 
i take 16 points at (3 = 6.0 as x = {x±, X2, x%, x±) with 
(a\ Xi S {0, l}{i — 1,2, 3), £4 € {0,nt}, and 2 points at 
(3 = 6.1 and 6.2 as x e {(0,0,0,0), {1,1,1, m)}- In ad- 
dition, the summation over the KS-hypercube is taken 
on each selected space-time point. At each quark mass 
m and temperature, we thus achieve high statistics as 
1600 data at f3 = 6.0 and 200 data at f3 = 6.1 and 6.2. 
Note that for the lattices of 20 3 x 8 (fi = 6.1) and 24 3 x 10 
(/3 = 6.2), we obtain 2000 data, which guarantees relia- 
bility of the results even in the vicinity of the critical 
temperature. 

At each temperature, we calculate the condensates (qq) 
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and g(qa ^ V G ' ^q) at three different quark masses, m. We 
observe that both the condensates show a clear linear 
behavior against the quark mass at every finite tempera- 
ture, as was seen in Ref. [l2] at T = 0. Therefore, we fit 
the data with a linear function and determine the con- 
densates in the chiral limit. 

We comment here on the error estimate and the numer- 
ical check on the finite- volume effect in lattice QCD. The 
statistical error is estimated with the jackknife method. 
In the vicinity of T c , i.e., on 20 3 x 8 at (3 = 6.1 and 
24 3 x 10 at ^ = 6.2, we conservatively suppress the auto- 
correlation by performing binning of the size of 10. We 
find that the statistical error is typically 5-7% level for 
every temperature. For the check of the finite- volume ar- 
tifact, we examine the condensates imposing the periodic 
condition on the KS-fields, \ and x, instead of the anti- 
periodic condition, at the 3-dimensional spatial bound- 
aries. We perform the calculations with typical lattices 
of 24 3 x 10 (J3 = 6.2) and 16 3 x 4 (/? = 6.0), which cor- 
responds to T ~ T c and T c -C T, respectively, and find 
that the dependence of the condensates on the spatial 
boundary condition is only 1% level. For T <C T c lat- 
tice case, it was shown that the finite- volume artifact for 
16 4 (/3 = 6.0) lattice is only 1% level in Ref.0. We there- 
fore conclude that the physical volume L 3 J> (1.6fm) 3 in 
our simulations is large enough to avoid the finite volume 
artifact. 



III. THE LATTICE QCD RESULTS AND 
DISCUSSIONS 

We estimate the thermal effects on each condensate, 
g(qa^ v G^q), or (qq), by taking the ratio between the 
values at finite and zero temperatures. In this ratio, the 
renormalization constants are canceled because there is 
no operator mixing for both the condensates in the chiral 
limit [17j. Here, at each j3, we use the value at the lowest 
temperature as a substitute for the value at zero tem- 
perature. This approximation can be justified because 
the finite-volume artifact for the lowest temperature is 
almost negligible, as was shown for 16 4 (/3 = 6.0) lattice. 

In figure^] we plot the thermal effects on g (qa ^ V G ^q) . 
We find a drastic change of g(qa^ v G^q) around 
the critical temperature T c ~ 280MeV. This is 
the first observation of chiral-symmetry restoration 
through g(qo ^ G^q). We obtain T c j ' = 0.64(4), 
which is consistent with the coincidence of the con- 
finement/deconfinement phase transition and chiral- 
symmetry phase restoration. We also find that the ther- 
mal effects on g(qa^ u G^ v q) are remarkably weak below 
T c , namely, T <, 0.9 T c . In figure^ the thermal effects on 
(qq) are plotted, and the same features as g(qa^ v G^ u q) 
are obtained. 

We then quantitatively compare the thermal effects of 
9{q\ a tivG fiuQ) and (qq). For this purpose, we plot the ratio 
m l( T ) = 9{^hvGhv1)t/ (<7<7>t normalized by mg(T = 0) 
for T < T c in figure [3] We observe that rn^(T) is almost 
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FIG. 1: The quark-gluon mixed condensate g(qa llv G^q)T 
at finite temperature T normalized by <7{g<7 M „G M „g)T=o- The 
vertical dashed line denotes the critical temperature T c ~ 
280MeV in quenched QCD. 
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FIG. 2: The quark condensate (qq)r at finite temperature T 
normalized by {qq) T=0 . The vertical dashed line denotes T c . 



independent of the temperature, even in the very vicinity 
of T c , which can be interpreted that giqc^G^qjx and 
(qq)T obey the same thermal behavior. This result is 
rather nontrivial because, as was noted before, these two 
condensates characterize different aspects of the QCD 
vacuum. 

2 2 

m Q (T) / m Q (T=0) 
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FIG. 3: The ratio rn^iT) = g(qa ^ V G ^q)T I (<jq)T normalized 
by mo(T = 0) plotted against T. This result indicates the 
same chiral behavior between g(q<J l _ lv G l _ lv q)T and {qq)r- 
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As a new interesting feature absent in the scalar-type 
condensate (qq), giqu^G^q) may reveal possible break- 
ing of the space-time duality in Euclidean QCD at fi- 
nite temperature. For further analysis in this direction, 
we investigate the thermal effects on electric/magnetic 
components of g(qo ^G p, v q) separately, i.e., g(qouGuq) 
and g(qajkGjkq) ■ Hereafter, no summation is taken for 
k € {1,2,3}. Of course, the equality g(q<JnG/nq) = 
g(qcrjkGjkq) holds at zero temperature from 0(4) sym- 
metry in Euclidean QCD. This equality, however, does 
not necessarily hold for finite temperature, where the 
space-time duality is explicitly broken. In fact, at finite 
temperature, g(qanGaq)T and g(q&jkGjkq}T can be re- 
garded as different chiral order parameters. 

We therefore plot the ratio Re/b(T) defined by 



Re/b(T) 



Y,g{qv4iGaq)T 

i 

E g{qO]kG Jk q) T 

j<k 



(5) 



at finite temperature T <J T c in figure One sees that 
the equality g{qcTuGuQ)T — g(q&jkGjkq)T holds for even 
at finite temperature. This again indicates that all the 
(/!, v) components of g(qo ^G ilu q)T obey the same chiral 
behavior near T r . 
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FIG. 4: The ratio Re/b{T) = 

'52 i g{qouGuq)T/Y J j <k 9 l \q (J 3kGjkq)T plotted against 
the temperature T. This result indicates g(qcr4,iGuq) and 
giqvjkGjkq) obey the same chiral behavior. 

In this way, we observe that the chiral condensates, i.e., 
(QQ)t, g(qcr^G fn ,q) T , g(q<J4iG 4 iq)T and g(qcr jk G jk q) T , 
obey the almost identical behavior near T c . One of the 
possible explanations of this similarity is given by the 
so-called "Swiss cheese picture". In this picture, phase 
structure of the QCD vacuum at the critical temper- 
ature is represented by the mixture of co-existing two 
phases, i.e., chiral-broken phase and chiral-symmetric 
phase. Since the deconfinement phase transition is of 
the first order in the quenched approximation, the co- 
existence of two phases is plausible. In this scenario, the 
common thermal behaviors of the chiral condensates can 
be explained if the thermal effects are mainly character- 
ized by the volume ratio of the two phases and the ther- 
mal effects on each phase are subdominant. Since this 



consideration can be applied to any chiral condensates, 
this scenario predicts a universal behavior of all chiral 
order parameters near T c , which provides new aspects on 
the chiral structure of the QCD phase transition. 

The common thermal behaviors of the chiral conden- 
sates can be also explained from a different viewpoint. 
Following Refs. [Til fli^ . we can express the condensates 
as 



(qq) 

g(qv^ G ^q) = 



= I [ d X> 



A' 2 



, mp(X') 



(6) 



where |A) denotes the eigenvector of the Dirac opera- 
tor as ip\X) = A|A), and p(A) the spectral density on 
A. When one takes the chiral limit m — > after the 
thermodynamic limit V — * 0, the zero modes (A = 0) 
of the Dirac operator are responsible for both the con- 
densates. If (A|<7 ^ V G \j, v \\) does not have a singularity at 

A = 0, TOq corresponds to (X\(J^G^ U \\) . The uni- 

A— 

versal chiral behavior against the temperature indicates 
that both the electric and the magnetic components of 

(A|ct M1/ G for the chiral zero modes have remark- 

A— 

ably weak dependence on the temperature for T < T c , 
even in the very vicinity of T c . 



IV. SUMMARY 

We have studied thermal effects on giqcr^G^q) , in 
comparison with (qq), using the SU(3) C lattice QCD 
with the KS-fermion at the quenched level. For T <, 
0.9 T c , both the condensates show very little thermal 
effect, while a clear signal of chiral restoration is ob- 
served near T c as a sharp decrease of both the conden- 
sates. It has been a surprise that the ratio m§(T) = 
g(q cr fj,vG l j J vq)T/ (qq)r, hi contrast to the drastic change of 
each, is almost independent of the temperature in the en- 
tire region of T up to T c . It is also found that the Lorentz 
components, g(qo4iGnq)T and g(q<jjkGjkq)T 1 show an 
identical temperature dependence. Such observations in- 
dicate that these chiral condensates show a common ther- 
mal behavior, which may reveal new aspects of the chiral 
restoration of the QCD vacuum. For further studies, a 
full QCD lattice calculation is in progress in order to 
analyze the dynamical quark effects on the condensates. 
Through the full QCD calculation and its comparison 
with quenched QCD, we expect to get deeper insight on 
the chiral structure of the QCD vacuum at finite temper- 
ature. 
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